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ABSTRACT 

The measurement of the forward-backward asymmetries of heavy quarks 
provides one of the most precise determinations of sin 2 6 W in Z decays. 
We discuss in detail the one-loop QCD radiative corrections to these 
asymmetries. Results are given for single heavy-quark jet asymmetries 
and asymmetries of the thrust axis, as well as for heavy-quark two-jet 
final states. 



1. Introduction 

The heavy quark sector, especially 5-quark physics, has received much attention in the high- 
precision analyses of the electroweak part of the Standard Model which are being performed 
in Z decays at LEP and SLC [for reviews see Ref. |]J]. In particular, the measurement of 
the forward-backward asymmetry for 6-quark production A b FB provides one of the most precise 
determinations of the electroweak mixing angle sin 2 9w m Z decays 0. The precision of the 
experimental data is so high that higher-order electroweak and QCD corrections must be 
incorporated properly 

The photonic and the genuine electroweak radiative corrections to quark forward-backward 
asymmetries Ap B have been discussed thoroughly in the literature, Ref. ||. Many features of 
the radiative corrections due to strong interactions have also been presented in the past [4- 
13]. At the one-loop level the QCD corrections consist of virtual quark-vertex corrections and 
bremsstrahlung corrections where one of the quarks emits an additional gluon in the final state. 
This bremsstrahlung correction will depend in a sensitive way on the quark-jet properties, i.e. 
the definition of the jet axis and the maximum invariant mass allowed for the final-state jets. 

In this paper we will give a systematic and complete overview of the one-loop QCD correc- 
tions to the forward-backward asymmetries for heavy-quark production in e + e~ annihilation 
near the Z resonance. We will complete the theoretical analysis and also correct some erroneous 
points in the literature. We will focus on the quark-mass effects; with increasing accuracy of 
the LEP experiments, these effects become more and more relevant for the production of bot- 
tom and charm quarks. 

The results will be presented for two different definitions of the reference direction: the heavy 
quark axis which is usually used in analytical calculations ["single quark-jet asymmetry"], and 
the thrust axis, corresponding to the direction of the most energetic jet, which is more relevant 
from the experimental viewpoint. We will also present results for forward-backward asymme- 
tries of two-jet final states in which the invariant mass of the jets is less than a fraction ^fy of 
the total energy. In these two-jet asymmetries, the QCD corrections are much smaller than for 
the total (single quark-jet and thrust-axis) asymmetries. Two-jet asymmetries are therefore 
very useful to minimize the errors due to strong interactions in high-statistics measurements. 

The paper is organized as follows. In the next section, we will summarize, for the sake of 
completeness, the tree-level results for the forward-backward asymmetries. In section 3, the 
QCD corrections will be discussed for the case in which the reference directions are either the 
quark or the thrust axes, and the results for the single quark-jet and thrust-axis asymmetries 
will be presented. In section 4, we will define the two-jet forward-backward asymmetries; 
results of the QCD corrections will be given for the experimentally useful range of invariant 
masses of the jets in the final state. A short summary in section 5 will conclude the paper. 



2. Basic set-up 

The differential cross section for the process e + e~ — > // is a binomial in cos#, where 6? is the 
angle between e~ and /, 

1 + cos 6')cr[/ + -sin tfo" L H — cos yap (1) 



dcos# 8 V y 4 4 

C/(L) denote the contributions of unpolarized (longitudinally polarized) j/Z gauge bosons along 
the reference axis, and F denotes the difference between right and left polarizations. The total 
cross section is the sum of ajj and o^, 

a T = (Tu + ctl (2) 

The forward-backward asymmetry for the full angular range can be expressed by the ratio of 

Op to Op, 

A' FB = \%- (3) 

4 (J? 



In the Born approximation, the parton cross sections <tt and op are given by 

a T = P^cr vv + p 3 a A A 

a F = p 2 a V A (4) 



with (5 = J\ — 4m 2 Js = y/1 — fi 2 being the velocity of the final quarks. [The quark mass, 
denoted by m/, will later be defined more precisely] The phase space suppression and the 
dynamical P-wave suppression are entirely taken account of by the velocity-dependent coeffi- 
cients. In terms of the electric charge of the fermion e/ and the vector and axial- vector couplings 
of the fermion to the Z boson [J| L = ±1/2 is the weak isospin, and s 2 ^ = 1 — c\y = sin 2 9w], 



v f = 2lf - As 2 w e f and a f = 2lf 
the cross sections, for quark final states, can be written as 

47ra 2 (s) 2 2 G F a(s) M 2 (s - M 2 ) G\ Mis , , 2x 2 

°AA = ^^T^ + a l) a ) ( 5 ) 

Gpa(s) M z {s - M 2 Z ) G 2 F M z s 

ctva = — 2 ^| e e e f a e a f + —^v e v f a e a f 

where, including the energy-dependent Z-boson width, D 2 Z is given by the Breit-Wigner form 
D\ = (s - M 2 ) 2 + (sT z /M z ) 2 . 

On top of the Z resonance, the forward-backward asymmetry is dominated by the pure Z 
exchange amplitude so that in the Born approximation 

f 3 2v e a e 2pv f a f 

FB/0 4 v 2 + a 2 e |(3 - P 2 )v 2 + P 2 a 2 { } 



3. Asymmetries of single quark jets and thrust axis 

At the one-loop level, the QCD corrections to the forward-backward asymmetry for quark- 
pair production can be divided into three parts: (a) the virtual vertex correction; (b) soft 
gluon bremsstrahlung [the additional gluon in the final state carries an energy less than A = 
2E g /^/s <C 1]; and (c) hard gluon bremsstrahlung [the gluon carries an energy larger than E g \. 
The results for these corrections in Z decays are summarized below. 

Virtual corrections and soft bremsstrahlung: 

The contribution of the virtual corrections to the forward-backward asymmetry, Ap B , is given 

by 

SA Q I __ A Q x 2a, W 2 a 2 Q ~H 1 - (3 2 1-/3 

^ FB | virtual - a fb/0 ^^l v 2 q + ^ q -p lo s TTp {7) 

The infrared singularities which are present in both ot and o>, cancel in the ratio so that 
5Ap B (virtual is free of any singularities. Note that for quark masses approaching zero, the 
virtual corrections vanish. 

The soft gluon bremsstrahlung in the process e + e~ — > QQg, in which the additional gluon 
carries an energy less than A = 2E g j ^fs <C 1, contributes to ot and op by the same amount, 

a s F f brems = a F , T (l + 5 SB ) (8) 

The analytical expression of <5sb can be found in Ref. for instance. This correction therefore 
cancels out from the ratio op/or, and soft gluon bremsstrahlung does not contribute to Ap B 
to first order in as, 

Vsip B \ so ft brems — ^FB/0 X \ ) 

The cancellation of the soft bremsstrahlung correction is physically plausible since the emis- 
sion of soft energy quanta cannot affect the angular direction of the leading heavy quarks, a 
consequence of Galilei's law of inertia. 

Hard bremsstrahlung: 

In a first step, when analyzing the semi-inclusive process 

e+e- -► QQg ^Q + X 

we assume the quark momentum to be measured precisely with the quark direction to be 
chosen as the reference axis. With this choice, the contribution of the hard bremsstrahlung to 
the forward-backward asymmetry in Z decays including mass effects, is given by 
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Here, x and x denote the quark and antiquark energies in units of the beam energy, z = 2 — x — x 
is the scaled gluon energy, x — Vx = (1 — x )/(l — ^) anc i ^± = X + X =t 2. p and p are the 
momenta of the quark and antiquark relative to the reference axis and p± is the component of 
the antiquark momentum perpendicular to this direction: 
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Using the more convenient variables z and w 

1/2 



X 



x the phase space integration has to be 



performed over \w\ < z [1 — ft /(l — z)\ and 2A < z < 1 — /jr. 

The first contribution in Eq. (13) accounts for the QCD corrections of the genuine vector- 
axial vector interference term while the second part is due to the corrections to the total cross 
section. 

However, due to fragmentation effects, cascade decays and detector imperfections, the quark 
direction cannot always be determined very well. In this realistic situation, the thrust axis lends 
itself to the reference axis of the forward-backward asymmetry. While the thrust axis coincides 
with the quark/antiquark line in quark-pair production to lowest order, this connection is 
destroyed by hard-gluon emission. In fact, with a small probability the gluon carries even more 
energy than the quarks/antiquarks in three-jet events. 

The experimental procedure of choosing the thrust axis and its orientation introduced in 
Ref. [|l4[], coincides with the definition of the reference axis in this theoretical analysis. The 
orientation of the thrust axis T is defined in such a way that T • p is positive where p is the 
3-momentum of the b quark (Fig. 1). 

In the three different regions of the Dalitz plot we have for the momenta p, p and p±: 

• x > x, z [reference axis parallel to quark line] 



P — \J X* — /l* , p = \jX 2 — /l 2 COS 9 q q , p± = \J X 2 — /J? SU1 9 q q 

x > x, z [reference axis anti-parallel to antiquark line] 



P 



x 2 — /i 2 



V 



X 2 — /l 2 COS 9gq , p± = 



z > x, x [reference axis anti-parallel to gluon line] 



p = \jx 2 - p 2 cos qg , p= \jx 2 - p 2 cos Ogg , p± = \j x 2 - p 2 sin 6q g 

with 

2(1 — x — z) + xz 2(1 — x — z) + xz 

cos 9qg = and cos 9 qg = . 

z^/x 2 — p 2 z^/x 2 — p 2 

In addition, Rp has to be multiplied by (—1) when the reference axis coincides with the an- 
tiquark or gluon axis [change from cos 9 to — cos 9}; in R v , the term xp\ has to be replaced 
by \+p 2 j_ when the gluon axis is chosen as the reference axis. Note that according to our pre- 
scription for the orientation of the thrust axis, the first two cases give identical results. Events 
in which the gluon 3-momentum coincides with the thrust axis, dilute the forward-backward 
asymmetry since / <jp integrated over the proper part of the Dalitz plot must vanish for this 
class of events. 

The sum of the virtual soft corrections and the hard gluon bremsstrahlung, integrated over 
the entire allowed range, will change the tree-level forward- backward asymmetry of single 
quark-jets to 



aQ __ aQ 

J± FB — -H-pB/Q 
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(11) 



The coefficients C F [which are free of any singularities and independent of the gluon energy 
cut separating the soft from the hard bremsstrahlung] are shown in Table 1 for the production 
of quarks with nearly zero masses, c-quarks and 6-quarks. Two reference directions have been 
chosen: the quark axis and the thrust axis. The numbers to the left of the bars correspond to 
the pole masses of the quarks, m c = 1.5 GeV and m^ = 4.5 GeV; the numbers to the right to 



the choice |T5| m c = 0.7 GeV and m& = 3 GeV for the c, b mass parameters, the mass values 
in the MS scheme at the scale Mz- The difference between the two predictions reflects the 
uncertainties due to the corrections of next-to-leading order which have not been included in 
this analysis. These uncertainties in the prediction of the FB asymmetries, of order 10~ 4 , are 
much smaller than the anticipated ultimate error of the LEP measurements. 

One first notices that the coefficient Cp for a given quark species is smaller when the 
reference axis is defined as the thrust axis. However, while the difference is about ~ 10% for 
light quarks it decreases to ~ 5% in the case of b quarks with mj = 4.5 GeV. The numerical 
results in the case where the quark axis is chosen as the reference axis, agree with those given 
in Refs. [5-7], and the result in the massless case for the thrust axis being the reference axis 
agrees with Ref. [110] if the pole masses are chosen. 

The prediction (— a s /n) for the QCD correction of the FB asymmetry in the massless case 
has been noticed first in Ref. ||. The simple coefficient (—1) is a consequence of the fact that 
the QCD corrections to ap vanish so that the FB asymmetry is only affected by the well-known 
QCD correction 1 + a s /n to the total cross section in the denominator. 

Note also that in the case where the quark axis is chosen as a reference, the coefficient Cp 
is ~ 20% smaller for 6-quarks with m& = 4.5 GeV than for massless quarks [as anticipated 



Ref. Axis 


/^m=0 


C F 


C F 


Quark 


1.00 


0.93 | 0.96 


0.80 | 0.86 


Thrust 


0.89 


0.86 | 0.88 


0.77 | 0.81 



Table 1: Coefficients Cp of (— a s /7r) in the single quark-jet forward-backward asymmetry with 
the reference axis being the quark axis, and for the thrust axis; left: pole quark masses, right: 
MS masses at the scale Mz- 

from Galilei's law of inertia]. Therefore, although //& = 2m&/Mz ~ 0.1, quark mass effects are 
vey important for the QCD corrections. For this choice of the reference axis, the coefficient Cp 
has first been calculated numerically in Ref. |5j for arbitrary quark masses, later it has been 
derived analytically || . A very good approximation can be obtained by making an expansion^ 
in the mass parameter /i || Fjj ; 
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In turn, if the thrust axis is chosen as the reference direction, the mass effects are slightly less 
pronounced. Note also that it is very difficult to obtain analytical results for Cp in this case 
due to complicated phase space integrals. 



4. Two— jet forward— backward asymmetry 

The two-jet asymmetry is defined for all events in which the invariant mass of the jets is less 
than a fraction sjy of the total energy, i.e. 

(p + p) 2 /s<y — 



(P + VgY/s < y 

(P + P 9 ) 2 /s <y 



z > 1 -y 
x > 1 +\fi 2 
X > 1 +i/i 2 



y 
y 



The two-jet forward-backward asymmetry may then be written again as 
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(13) 



x Due to a small error in the analytic expressions of the integrated parton cross section of Ref. pl| , the 
approximate formula for the forward-backward asymmetry obtained in Ref. JTJ] was slightly incorrect. The 
result was so close to the correct prediction that numerical cross-checks failed to reveal this tiny discrepancy. 
The effect on the asymmetry was of order 10 -4 for b quarks and thus about two orders of magnitude below the 
present experimental error, for c quarks smaller still by another order of magnitude. 



with the coefficients Cp given in the Table 2 and 3 for a sample of y values^]; the reference axis 
is chosen to be the quark axis in Table 2 and the thrust axis in Table 3. The quark masses 
have been defined as in Table 1. The entries to the left of the bars again correspond to the pole 
quark masses, those to the right to the MS mass definitions at the scale Mz- 



Z/cut 


f~tm=0 


c c F 


cp 


0.01 
0.02 
0.04 
0.08 
0.16 


0.12 
0.21 
0.35 
0.55 
0.81 


0.06 | 0.09 
0.14 0.17 
0.28 0.31 
0.48 0.52 
0.74 0.77 


0.02 | 0.03 
0.05 0.08 
0.14 0.20 
0.34 0.40 
0.60 0.66 


2/3 


1.00 


0.93 | 0.96 


0.80 | 0.86 



Table 2: Coefficient Cp of (—o> s /ir) in Ap B as a function of the invariant mass of two jets; the 
reference axis is taken to be the quark axis; left: pole masses, right: MS masses at the scales 
M z . 



Z/cut 


/~tm=0 


C F 


cp 


0.01 
0.02 
0.04 
0.08 
0.16 


0.07 
0.13 
0.23 
0.39 
0.64 


0.04 | 0.06 
0.10 0.11 
0.19 0.21 
0.35 0.37 
0.60 0.62 


0.02 | 0.03 
0.04 0.06 
0.11 0.15 
0.26 0.30 
0.51 0.55 


2/3 


0.89 


0.86 | 0.88 


0.77 | 0.81 



Table 3: Coefficient Cp of (—a s /ir) in Ap B as a function of the invariant mass of two jets; the 
reference axis is taken to be the thrust axis; left: pole masses, right: MS masses at the scale 
M z . 



The following conclusions can be drawn from the tables. First, as anticipated, the two-jet 
coefficient Cp decreases with decreasing y cut, and they become very small for small invariant 
masses. Two-jet asymmetries are therefore less affected by the QCD corrections and they 
can be used to reduce the errors coming from the relatively poor knowledge of a s and from 



2 The entries in Tables 2/3 supersede those given in Ref. |7J]. 



higher-order effects. Furthermore, the coefficient Cp of the two-jet asymmetry follows the 
same pattern as the asymmetry for single quark jets and for the thrust axis in Table 1: The 
QCD corrections are smaller in the case where the thrust axis is chosen as the reference axis, 
than in the case where the reference axis is the quark axis; the corrections decrease also with 
increasing mass. Finally, for the maximum value of the y cut, y = 2/3, there are by definition 
no three-jet events left over and the two-jet asymmetry approaches the single-jet forward- 
backward asymmetry. 

5. Summary 

We have presented an overview of the one-loop QCD corrections to forward-backward asym- 
metries for heavy-quark production in Z decays. A complete set of results has been given 
which improve the theoretical predictions for these asymmetries, the measurement of which 
provides some of the most precise determinations of the electroweak mixing angle sin 2 9w in 
Z decays at LEP and SLC We have focussed on the consequences of quark-mass corrections 
which we have shown to be rather important. The difference between the QCD corrections to 
the forward-backward asymmetries for light quarks and for 6-quarks turns out to be as large 
as 20%. With the high-accuracy of the measurements performed at LEP, these effects must be 
taken into account properly. 

We have discussed these radiative corrections for the two cases where the reference directions 
are taken either to be the quark axis or the thrust axis. In the second case the QCD corrections 
are slighty smaller. We have also considered the consequences of constraints on the jets by 
demanding their invariant mass to be less than a fraction <Jy of the total energy. These two-jet 
asymmetries are affected less by QCD corrections than the single-jet asymmetries, and they 
could therefore be exploited to minimize the errors due to strong interactions in experimental 
high-statistics analyses. 
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Figure Caption 

Fig. 1 Definition of the oriented thurst axis T in 3-jet events. 
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